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ABSTRACT 

We use the spatially-resolved, multi-band photometry in the GOODS South field acquired by the 
CANDELS project to constrain the nature of candidate Lyman continuum (LyC) emitters at redshift 
z ^ 3.7 identified using ultra deep imaging below the Lyman limit (1-sigma limit of «30 AB in a 
2" diameter aperture). In 18 candidates, out of a sample of 19 with flux detected at >3-sigma level, the 
light centroid of the candidate LyC emission is offset from that of the LBG by up to 1.5". We fit the 
SED of the LyC candidates to spectral population synthesis models to measure photometric redshifts 
and the stellar population parameters. We also discuss the differences in the UV colors between the 
LBG and the LyC candidates, and how to estimate the escape fraction of ionizing radiation (/ esc ) in 
cases, like in most of our galaxies, where the LyC emission is spatially offset from the host galaxy. 
In all but one case we conclude that the candidate LyC emission is most likely due to lower redshift 
interlopers. Based on these findings, we argue that the majority of similar measurements reported 
in the literature need further investigation before it can be firmly concluded that LyC emission is 
detected. Our only surviving LyC candidate is a LBG at z = 3.795, which shows the bluest (B — V) 
color among LBGs at similar redshift, a stellar mass of M ~ 2 x 10 9 M Q , weak interstellar absorption 
lines and a flat UV spectral slope with no Lya in emission. We estimate its f esc to be in the range 
25%-100%, depending on the dust and intergalactic attenuation. 

Subject headings: galaxies: distances and redshifts - galaxies: high-redshift - galaxies: formation 



1. INTRODUCTION 

The origin of the ionizing radiation, whether from 
galaxies or AGN, responsible for the re-ionization of the 
universe at high redshift, z > 6, and for keeping it ionized 
at later epochs, is still poorly constrained. The contri- 
bution of quasars to the hydrogen ionizing background 
increases as we look back in time from z = to z ~ 2 as 
the peak of the quasar luminosity function is approached 
(e.g., Fanidakis et al. 2011). Beyond redshift 2 their con- 
tribution significantly decreases and the hydrogen pho- 
toionization rate is most probably maintained by stellar 
emission (e.g., Siana et al. 2008; Faucher-Giguere et al. 
2009; Haardt & Madau 2011; but see Fiore et al. 2011). 
The star-forming galaxies at z > 3 are therefore the lead- 
ing candidates to provide the remaining ionizing photons. 
The severe IGM absorption at z > 4.5 prevents us from 
observing directly the Lyman continuum (LyC) emission 
(Meiksin 2006; Inoue & Iwata 2008). This is particularly 
true during the re-ionization of the Universe (z > 7). It 
is therefore essential to identify diagnostics of the LyC 
emitters at lower redshift, z ~ 3 — 4 (at \ re st > 1216A), 
and infer if galaxies with these characteristics are more 
common during the epoch of re-ionization (one Gyr ear- 
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lier, or z > 7). 

At z ~ 3 — 4 the direct detection of LyC emission 
from galaxies is still difficult, because ionizing radiation 
is severely attenuated by neutral gas and by dust in the 
interstellar and circumgalactic medium (ISM and CGM) 
of the source itself, as well as by the intervening IGM. 
As a result, the number of reported detections of star- 
forming galaxies at high redshift, typically selected as 
Lyman-break galaxies or with equivalent criteria, with 
LyC emission is very small (Shapley et al. 2006; Iwata 
et al. 2009; Nestor et al. 2011; S06, 109 and Nil, here- 
after; Vanzella et al. 2010b, VlOb, hereafter; Boutsia 
et al. 2011). The issue is further complicated by the 
relatively high probability of finding faint low-redshift 
galaxies at very close angular separation, i.e., ~ 1 arcsec 
or less, from a high-redshift galaxy with brighter appar- 
ent magnitude (Vanzella et al. 2010a, VlOa, hereafter). 
These interlopers, for which spectroscopic identification 
is not viable, can be wrongly interpreted as spots of LyC 
emission from the host galaxy if high-angular resolution 
multi-band photometry is not available, as is often the 
case, to recognize their nature. 

Recently, searches of LyC emitters with deep imag- 
ing below the Lyman limit have yielded samples of can- 
didates at redshift z^ 3 where the region of emission 
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of the candidate LyC ionization radiation is spatially 
offset from that of the non-ionizing rest far-UV light 
(dubbed "nLyC" in what follows), namely the main 
body of the galaxy typically observed between 1500 A and 
2000A. The displacement of these putative LyC-emitting 
"blobs" is generally less than one arcsec, but in some 
cases values as high as 2" have been reported (109; Nil), 
corresponding to separation in the range of several kpc 
to a few tens of kpc. In other words, the candidate LyC 
emission would come from regions of the galaxies that are 
well separated from the center, which is where most of the 
LyC photons are produced, and in some cases so far from 
it to be classified as separate sources. If confirmed this 
would provide important empirical constraints on how 
LyC escapes from galaxies. Future observations with the 
new HST/WFC3 in the LyC rest-frame at z ~ 3 will 
help exploring this issue further. In the meantime, how- 
ever, motivated by the potential importance if the above 
findings are confirmed, we use spatially resolved, multi- 
band photometry from HST to constrain the nature of a 
sample of LyC emitter candidates that we have selected 
in the GOODS South, one of the fields targeted by the 
CANDELS project (Grogin et al. 2011; Koekemoer et al. 
2011) from ultra-deep U-band imaging below the Lyman 
Limit of star-forming galaxies with known spectroscopic 
rcdshift z > 3.4. We also report on the only one Ly- 
man Break Galaxy in our sample at z = 3.795 (Ionl, 
hereafter) with a clear LyC emission. 

2. ANALYSIS OF A SAMPLE OF CANDIDATES LYC 
EMITTERS IN GOODS-SOUTH 

We have selected a sample of LyC emitter candidates 
from ultra-deep imaging in the U-band of the GOODS- 
South field (VlOb). The U-band images, which are de- 
scribed in Nonino et al. (2009), were aimed to probe the 
spectral region below the Lyman limit of these galaxies 
for possible emission of ionizing radiation. The images 
are very deep, reaching 1-a flux upper limit of about 
30 mag (AB) for an unresolved source within a circular 
aperture of 2 diameter. 

A relatively large number of sources, selected as 
Lyman-break galaxies, have secure redshift identifica- 
tion in the GOODS-S field (Vanzella et al. 2008, 2009; 
Popesso et al. 2009; Balestra et al. 2010). As described 
in VlOb, we start from a spectroscopic sample of 122 B- 
band dropouts at 3.4 < z < 4.5 selected by Giavalisco 
et al. (2004) from the GOODS/ ACS images for which 
we have robust redshifts, i.e., Quality Flag QF=A. The 
lower limit of the redshift range is the lowest value such 
that the 912A Lyman limit is outside, redward of the 
system throughput in the U filter, while the upper is 
chosen because the rapid increase of opacity of the IGM 
produces too small transparency at higher redshift to 
make analysis of z > 4.5 galaxies useful (see VlOb). 7 
out of 122 B-band dropouts are detected in the Chandra 
4 Msec images in the CDF-S (Xue et al. 2011) and clas- 
sified as AGN. Among these 122 sources we have flagged 
as potential LyC emitter candidates the 32 of them which 
have flux detection at the 2 a level or larger within a cir- 
cular aperture of 1.2" diameter in the U-band image. 
In this work we discuss in detail these U-band emitters 
(named Vera hereafter). 

As we noticed in VlOb, in 28 of the 32 sources the U- 
band emission, i.e., the candidate LyC light, is spatially 



offset relative to the centroid of the rest-frame far-UV 
light of the LBG (the nLyC at wavelength around « 1500 
A). A visual inspection of the HST/ ACS images showed 
that in 9 of the 28 sources the U-band flux light can 
unambiguously be explained as coming from the outer 
isophotes of nearby sources, which most likely are fore- 
ground interlopers (e.g., see Figures 4 and 5 of VlOb). 
We have eliminated these galaxies from the sample of 
LyC candidates. The remaining 19 candidates are listed 
in Tab. 1 and are analyzed here. They have counterparts 
in the ACS and WFC3 images which are close to the tar- 
geted host LBG with angular separations in the range 

0. 4"-1.9". Image cutouts for these sources are shown 
in VlOb, here we show in Figure [2] and the Appendix 
the more critical cases where the angular separation is 
smaller than 1". All but one of the U-band sources 
are fainter, in AB magnitudes, than their ACS optical 
(rest UV) nLyC counterparts. Tab. 1 reports the ratio 
of f 1500/ f Ly c for the 19 objects. 

The remaining 4 U-band detections of the initial sam- 
ple of 32 are co-spatial with the ACS and WFC3 images 
of the host sources, in the sense that the centroid of the 
U-band light falls within the errors (~ 0.1") on the lo- 
cation derived from the optical HST/ACS images. As 
discussed in VlOb (and reported in their Table 3), 3 out 
of 4 are detected in the 4 Msec X-ray Chandra images 
(0.5-8 keV), and also show typical signatures of AGN 
in their optical spectra, such as high ionization emission 
lines (e.g., C IV, N V). The forth source, which we call 
Ionl, is our only robust candidate stellar LyC emitter 
and we will describe it further in Sect. 6. 

2.1. The Effect Of The Intergalactic Attenuation 

The method we are adopting here is based on an 
intermediate-band filter (FWHM = 350A, U-band), 
and is similar to the typical narrow-band imaging only 
in the cases where 3.4 < z < 3.5 (9 out of 19 LBGs), 

1. e. the AA between the Lyman limit of the galaxy and 
the red cut-off of the filter is minimal. In this case, the 
greater depth of our image compensate for the larger 
noise due to the broader bandpass. For example, if 
compared with the NB3640 narrow-band filter used by 
Nil (FWHM = 100A), the depth of our U-band image 
nearly exactly compensate for the differet filter setup. 

Differences between the methods arise if we consider, 
as in our approach, a variable redshift. Indeed, the LyC 
region probed by our filter is redshift dependent, e.g., it 
is A rest < 909,800,74lA at z = 3.4,4.0,4.4, respectively. 
We have investigated the effect of the IGM attenuation as 
a function of redshift by running MC simulations as in 
VlOb (where the IGM attenuation is extensively taken 
into account). To this end, as a reference, we adopt 
a magnitude of 2775 — 24 and a ratio /1500//L1/C = 7 
(VlOb, Siana et al. 2007). We then apply 10000 differ- 
ent IGM transmissions convolved with the U-band filter 
and add photometric noise to the estimated flux in the 
U-band (see (VlOb)). In these conditions, at z = 3.4 we 
retrieve 82l^% of the sources in their LyC (U-band) at 

S/N > 2 (88±l% if 2775 = 22). As redshift increases, 
the fraction of recovered sources decreases because the 
IGM attenuation (see Figure [J) . The inner box of Fig- 
ure [T] shows how the IGM affects the fraction of recov- 
ered sources as a function of redshift (normalized to the 
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z = 3.4 case). 

It is worth noting that the majority of the sources an- 
alyzed here are at z < 3.8 (15 out of 19). Moreover, two 
examples at relatively high redshift have recently been 
reported: a LBG at z = 3.8 (lonl, also discussed in this 
work) and a AGN at z = 4.0 with 2775 = 26.09, for which 
a LyC emission at A res t < 830A and A res t < 800A is de- 
tected at S/N of 5.2 and 3.3 by our deep U-band imag- 
ing, respectively (VlOb). This suggest that despite the 
low average IGM transmission at z > 3.8, its stochas- 
tic behavior makes the variance to be relatively large. 
Indeed, at fixed redshift, the distribution of the IGM 
transmissions (U-band convolved) is asymmetric with 
an extended tail toward high values (Inoue & Iwata 
2008). The reason is that the Lyman continuum ab- 
sorption by the IGM is very stochastic because it is re- 
lated to the presence of relatively rare Lyman limit sys- 
tems (LLS) or damped Lyman a system (DLA), having 
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transmission is suddenly cut down at the corresponding 
wavelength. Conversely, without a LLS (or DLA system) 
near the source, we can expect a significant transmission 
even far below the source Lyman limit (see Inoue & Iwata 
2008). 

Having this in mind, we decided to keep the whole 
sample up to z ~ 4.4 in the analysis performed in the 
present work. If we could establish that the U-band 
selected sources observed in proximity of their LBG were 
at the same redshift of their companions then we would 
safely conclude that they are LyC emitters. Before to 
investigate the nature of their redshift, in the following 
section we briefly recall the issue of the contamination by 
lower redshift sources randomly placed at small angular 
separations from the higher redshift LBG. 

3. THE OCCURRENCE OF FOREGROUND 
CONTAMINATION 

The likelihood that a foreground interloper located in 
the vicinity of a LBG is responsible for the detection in 
the U band increases with redshift of the LBG, with the 
angular separation and the quality and depth of the ob- 
servations. Siana et al. (2007), 109 and Nil calculated 
analytically the probability that candidate LyC emission 
from a LBG is due to an interloper, given the character- 
istics of the data (quality and depth). VlOa performed 
the same calculations and ran Monte-Carlo simulations 
to quantify this effect. 

The displacement of the U-band light relative to that 
at redder wavelengths, as measured from the light cen- 
troid in the U-band and ACS images, for the sample 19 
galaxies discussed here, is in the range 0.4 < A9 < 1.9 
arcsec; in all cases there always is a counterpart to the 
U-band source in the HST/ ACS and WFC3 images. We 
calculated that the number of sources observed in the two 
annular bins with radii 0"-l" and l"-2" from the LBG 
centroid in the ACS 085O _ band images is equal, within 
the errors, to the expectations for foreground galaxies at 
increasing separations. That is, given the number counts 
and assuming a uniform distribution, the fraction of in- 
tercepted foreground sources increases with the area of 
the annulus considered. 

One of the parameters adopted in VlOa was the seeing 
(the PSF of the images), strictly related to the possibility 
to deblend close sources and related to the probability of 



superposition, which increases for the worst seeing con- 
ditions. The fact that the LyC emission could be intrin- 
sically offset from the main galaxy (as reported in Nil 
and 109) further complicates the interpretation, if the 
redshift of these potential emitters is not known. This 
is still true if high spatial resolution images are available 
and the sources are well separated, i.e., in absence of the 
redshift information an intrinsically offset LyC emission 
is fully compatible with the emission of a lower redshift 
object. 

On the one hand, the effects of this foreground con- 
tamination can be corrected statistically. On the other 
hand, it is worth investigating carefully each LyC candi- 
date, since any consideration about the mechanisms that 
allow the LyC photons to escape primarily depend on the 
reliability of the LyC detection. Therefore we now turn 
to the discussion of observational evidence that will help 
us to constrain the redshift, and hence the nature, of 
these sources in our sample (Sect. 4). 

4. CONSTRAINING THE REDSHIFT OF THE CANDIDATES 
LYC EMITTERS IN THE GOODS-S SAMPLE 

4.1. The Escape Fractions 

In this section we calculate the absolute and relative 
Lyman continuum escape fractions (defined below) of our 
potential LyC emitters assuming that they are at the 
same redshift of the LBG. Since f esc and (f esc ,rel) have 
to obey to clear limits, this in turn puts constraints on 
other quantities, in particular the redshift. In order to 
do that, we have to first define the relation among var- 
ious quantities and calculate the escape fraction where 
spatially offset LyC and nLyC emissions are present. 

4.1.1. Escape Fraction From a Morphologically Resolved 
Lyman Continuum 

Following Siana et al. (2007) the observed flux ratio 
between the 1500 A and the Lyman continuum is affected 
by several factors and is expressed as : 
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where LyC is the wavelength at which the Lyman con- 
tinuum is observed, (L1500/Ll v c)int is the intrinsic lu- 
minosity density ratio, (Ai 50 o — ^L y c) is the differential 
dust attenuation (in magnitudes), t hi.ig m {LyC) is the 
line-of-sight opacity of the IGM for LyC photons (and 

transmission is Tj^^f — e~ TL y°), and THi,iSM(LyC) 
is the optical depth of the Lyman continuum absorp- 
tion from H 1 within the observed galaxy's interstellar 

pHI _ 
ISM — 



medium, ISM (whose transmission is defined as T^ 1 



-THi,iSM(LyC)> 

The relative fraction of escaping LyC photons relative 
to the fraction of escaping nLyC (1500 A) photons is ob- 
tained rearranging the above equation : 



_ (L1500/L LyC )int , IGM-, 

/cscrci- ( ^i500/^ c ) obs eXP ^ cJ 



(2) 



it compares the observed flux density ratio (corrected for 
the IGM opacity) with models of the ultraviolet spectral 
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energy distribution of star-forming galaxies. If the dust 
attenuation Ai$oo is known, / OS c,rei can be converted to 
/osc as .f osc = lCr - 4Al500 / CSC:rc i (e.g., Inoue et al. 2005; 
Siana et al. 2007). Again, from the above equations, J esc 
can be written as: 

jcsc = exp [-T HI , ISM (LyC)] x i -°- 4 (^^), (3) 

the two factors on the right side have values in the range 
[0-1]. Clearly, their product, i.e., / CSC j cannot be greater 
than 1. 

It has been recently argued by Nil that due to the 
details of the radiative transfer and sources of the corre- 
sponding photons, any Lya emission and escaping LyC 
flux will not be necessarily co-spatial with either each 
other or with the bulk of the rest-frame UV flux in a given 
galaxy. This is true when considering Lya emission, 
which typically arises from backscattering of moving hy- 
drogen gas that can be spatially decoupled from the ion- 
izing sources. The observed light centroids (barycenter) 
of LyC and nLyC can be displaced in the case in which 
the LyC arises from a sub-region of a larger area, how- 
ever, the local emission is not spatially shifted, i.e., if the 
ionization radiation is measured in some sub-region, then 
the nLyC radiation is expected to be detected too (typ- 
ically with brighter magnitude). The crucial point here 
is that the two quantities (i*T500) o b s md (Fh y c)ohs must 
be measured in the same spatial (i.e., physical) region, 
where the ionizing and non-ionizing radiation arise. 

Conversely, recent works have performed measures 
with the aim to include both the fluxes (F1500) o b s and 
{FLyc)obs by enlarging the apertures (109) or by deriving 
total fluxes from the images regardless of the misalign- 
ment (e.g., SExtractor MAG- AUTO in U and R bands, 
as in Nil), resulting in measures within different aper- 
tures size and/or shape. In this way, the measured ob- 
served ratio (-F1500/Fij / c , )obs, and consequently the / esc 
quantity, are strongly biased. In particular the resulting 
/esc would be severely underestimated, because the cor- 
rect (F1500) o bs value would be smaller. Figure [3] shows 
an illustrative example in the GOODS-S field (also dis- 
cussed below and reported in Tab. 1) in which the ob- 
served ratio has been calculated with and without taking 
into account of the offset emission, top and bottom pan- 
els, respectively. Under the same assumptions (intrin- 
sic luminosity ratio and IGM transmission), the former 
method (top panel) produces an / csc that is ~ 10 times 
lower than the latter (bottom panel). 

The consequence is that if the estimated / esc exceeds 
the value of 100% then some other quantity has to be 
revised. It is even more significant if this happens under 
conservative assumptions of (L1500/Li y c)int, dust and 
IGM attenuation. 

If, from one hand, the / osc has to be less than 100%, 
the constraints can be even more stringent if / esc ,rei has 
to be less than 100%. Indeed from Eq. [T]and [5] it turns 
out that : 

1Q -0A(A 15OO ~A LyC ) x f = exp[-T HIJS M(LyC)} 

(4) 

In order to have the transmission of the ISM correctly in 
the range [0-1] (right side), the f esc ,rei has to be less than 
IqOA(a 1500 -a Lv c) ; w hich in turns is a quantity less than 



one (see Siana et al. 2007, their Figure 2). Therefore 
also the relative escape fraction has to be less than one. 

4.1.2. Anomalous Escape Fractions In Our Sample 

Following the discussion of the previous section, we 
have estimated / CS c,rei for our sample in the correct way, 
i.e., in same spatial regions where the [/-band emission is 
observed, assuming that the U-band sources are at the 
same redshift as the LBG, i.e., they are LyC emitters. 
Following Eq. [2] we calculate the / CS c,rci by assuming an 
intrinsic ratio (L1500/LL y c)int — 3 (Shapley et al. 2006) 
and the maximum IGM transmission at the redshift of 
the LBGQ The intrinsic ratio is also justified by the fact 
that the sources show an UV spectral slope (from their 
(i — z) color) similar or redder than their LBGs. There- 
fore the adopted ratio is a conservative assumption. 

The /esc,rci values are reported in the last column of Ta- 
ble 1. All but one U-band LyC candidates have / eS c,rei 
larger than 100%, that would implies a less dust atten- 
uation for shorter wavelengths, i.e., Aj, y c < A n L y c, in 
contrast to any dust extinction law. This is even more 
evident if we adopt the ratio (L1500/LL y c)int — 7 (e.g., 
Siana et al. 2007; VlOb), that increases the / eS c,rei 
values of Table 1 of a factor 7/3. The one case with 
/csc.rci < 100% has inconsistent photometric redshift and 
UV colors as described below. 

4.2. The UV Colors 

The B — V color of galaxies at the redshift considered 
here largely depends on the cosmic opacity of the IGM 
(see Madau 1995), specifically due to the effects of the 
Lya forest. Given the spatial correlation scale of the 
IGM, if the U-band companions are at the same red- 
shift as their associated LBGs, their B — V colors should 
therefore be similar. 

Figure H] shows the difference of the observed (B — V) 
color A(B — V) between the LBG and their compan- 
ions LyC candidates (empty circles, A(B — V) = (B — 
V)lbg~{B — V)uem ) as a function of the transverse sep- 
aration (physical) calculated at the mean redshift of the 
Lya forest. The same is shown for the four 4.0 < z < 4.5 
galaxies for which the A(V — i) has been adopted, more 
suitable than the (B — V) color in probing the IGM decre- 
ment. For all, the typical error in the color is <0.1 mag, 
because these sources are detected at more than 10-sigma 
in each band. As the figure illustrates, the LyC candi- 
dates are much bluer than their companion LBGs, im- 
plying either a very high fraction of escaping ionizing ra- 
diation and/or an high transmission of the IGM, or more 
likely as we are about to argue, that their redshift is sig- 
nificantly lower than that of the LBG, in agreement with 
the previous result and the photometric redshift analy- 
sis (reported in the next section). This is seen through 
a Monte Carlo simulations in which we have simulated 
f esc — 1 for the LyC candidates in the UV rest-frame 
continuum probed by the i?435-band and we have conser- 
vatively assumed that the IGM transmission at A < 912 

1 ^LyC 1 ^ as been calculated following the Inoue & Iwata (2008) 
prescription convolved with the VIMOS/U band filter, as described 
in VlOb; here we include in the calculations the recent statistics 
of the LLS provided by Prochaska et al. (2010) and Songaila & 
Cowie (2010), that increase slightly (see Inoue ct al. 2011). 

Transmissions at the redshift of the LBG have been calculated on 
10000 random line of sights. 
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A is the same as that at 912 < A < 1215 A (for the sim- 
ulations we used the IGM transmission by Inoue et al. 
(2008, 2011) and the SB99 UV templates with age 100 
Myr, continuous star formation and metallicity z=0.004 
and Salpeter IMF). We can then correct the observed 
(B — V) color of the LyC emitters and recalculate the 
A(B — V), which we plot in the Figure 2]as filled circles 
(case in which the f esc = 1). The corrected color differ- 
ences still show that the LyC candidates are much bluer 
than their putative LBG companions. An even much 
more conservative correction can be made by assuming 
unitary IGM transmission blueward of the Lya line, in 
all cases and recomputing the color differences (filled tri- 
angles), which still shows that the candidates would be 
much bluer than the LBG. 

The observed median (B — V) color and its 68% per- 
centile interval of the LyC candidates is 0.24^0 22) S1 &~ 
nificantly bluer than that of the LBG sample, which is 
1.85^0 5g, in contrast with what would be expected if 
their redshift were the same as that of the LBG and their 
stellar populations were similar (e.g., Meiksin 2006). 

We also explored a more extreme possibility by cal- 
culating at a given redshift and f esc — 1 the expected 
colors arising from population III stars (Pop III). The 
Pop III stars of Inoue (2011) template (with zero metal- 
licity), convolved with the U-band filter and the IGM 
transmissions (1000 lines of sight) are partially able to 
reproduce the observed (B — V) colors only in the case 
of high transparency of the IGM. Moreover, the presence 
of Pop III stars produces a very steep UV spectral slope 
(with a color (i— z) = —0.2), that is in contrast to what 
observed in our sample, i.e., < (i — z) > = 0.1 ± 0.2 0. 
Therefore we exclude the presence of Pop III stars. 

Regardless of theoretical expectations, in all cases that 
we considered the U-band companions have (B — V) 
color that are bluer by more than one magnitude (see 
A(B — V) in Tab. 1) than their corresponding LBG. If 
placed at the same redshift, the line of sight to each LyC 
candidate/LBG pair would be probing the IGM at trans- 
verse separations smaller than 20 kpc (physical) at the 
mean redshift of the Lyman forest. This is much smaller 
than the transverse cross-correlation function observed 
in QSO pairs separated of several arcminutes, for which 
the coherence length in the IGM is at least of the or- 
der of 500/i _1 kpc proper, (e.g., Fang et al. (1996); 
D'Odorico et al. (1998); Rauch et al. (2005); Cappetta 
et al. (2010)). This implies that the IGM attenuation 
due to the forest should be the same for the LyC candi- 
dates and their associated LBG if they were at the same 
redshift. Since the (i — z) colors of the Uem sources are 
similar or even redder than those of the LBG, imply- 
ing comparable or even redder ultraviolet spectral slope 
(similar stellar populations and obscuration properties), 
their (B — V) colors should also be similar, since the IGM 
attenuation would be in practise the same. 

Thus, we interpret these simulations as evidence that 
the LyC candidates are not regions of the LBG where 
ionizing radiation is escaping from but rather relatively 
unobscured star-forming galaxies at significantly lower 
redshift than the LBGs. 

2 The bluest sources in our sample (two with —0.2 < (i — z) < 
—0.1) have a (B — V) color that is bluer than the bluest value 
among the 1000 realizations calculated for Pop III stars. 



4.3. Photometric Redshift and SED Fitting 

To further investigate the nature of the U-band 
sources, namely whether regions of the host LBG from 
where LyC ionizing radiation is escaping or sources lo- 
cated at different, very likely lower, redshifts, we have 
fit the spatially-resolved CANDELS multi-band HST 
(BVizJYH), ground-based (U and K) and Spitzer/IRAC 
photometry to templates and spectral population syn- 
thesis models to derive their photometric redshifts and 
the parameters of their stellar populations (stellar mass, 
star-formation rate, dust obscuration and age) . 

The images in the above bands have different angu- 
lar resolution, and to robustly measure the photometry 
of our sources we have used the TFIT software package, 
developed by the GOODS and CANDELS teams (Lai- 
dler et al. 2007), which uses positional priors and PSF 
information to measure apparent magnitudes in matched 
apertures. For each source, TFIT uses the spatial posi- 
tion and morphology from the HST high-resolution im- 
ages (we used both the ACS z-band and the WFC3/IR 
H band) to construct a template image, which includes 
nearby sources. While these are fully resolved at the 
HST angular resolution, they might be marginally or 
even heavily blended in the other images. This template 
is then fit to the images of the object in all other low- 
resolution bands after convolution with the appropriate 
PSF. During the fitting procedure, the fluxes of the ob- 
ject in low-resolution bands are left as free parameters. 
The best-fit fluxes are considered as the fluxes of the ob- 
ject in low-resolution bands. These procedures can be 
simultaneously done for several objects which are close 
enough to each other in the sky. Experiments on both 
simulated and real images show that TFIT is able to 
measure accurate isophotal photometry of objects to the 
limiting sensitivity of the image (Laidler et al. 2007). 

We derive the photometric redshift and the physical 
properties of the LyC emitter candidates by fitting the 
observed spectral energy distributions (SEDs) to stellar 
population synthesis models. Models used to measure 
photo-zs are extracted from the library of PEGASE 2.0 
(Fioc & Rocca-Volmerange 1997). Instead of using the 
redshift with the minimum x 2 , we integrate the proba- 
bility distribution function of redshift (zPDF) and de- 
rive the likelihood- weighted average redshift. When the 
zPDF has two or more peaks, we only integrate the main 
peak that has the largest power. 

Since the U-band selected LyC candidates have coun- 
terparts in the HSTz and H bands in each case, we have 
generated two sets of fits, one based on positional priors 
from the ACS z-band images, which sample the rest far- 
UV SED of the galaxies at A < 1800 A, and the other on 
priors from the WFC3 H-band images (A < 3600 A). The 
first (GOODS) set uses photometry in the U and the ACS 
BViz bands, in the VLT/ISAAC JHK images of GOODS- 
South described by Giavalisco et al. (2004) and Rctzlaff 
et al. (2010), and in the four GOODS Spitzer/IRAC 
bands. The second (CANDELS) set uses the WFC3/IR 
YJH ncar-IR photometry in place of the ground-based 
one. The z-band positional priors might better reflect 
the spatial location of the LyC emitting regions, since es- 
sentially the same stars that power it also power most of 
the light in the nLyC far UV spectrum. The z-band im- 
ages also have better angular resolution than the H-band 
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ones. The second set offers an independent set of mea- 
sures, which takes advantage of the generally more sensi- 
tive and higher-resolution WFC3 data; it lacks, however, 
the K-band data. 

Since the potential emission in the LyC spectral range 
is not taken into account in the models, and thus flux 
in the U-band could skew the photometric redshift cal- 
culation towards lower values, for both photometric sets 
we have run the fit with and without the U-band pho- 
tometry. Figure [5] shows the redshift probability func- 
tion derived from the GOODS photometric set; Figure |6] 
the one from the CANDELS set. In all cases the solid 
curve refers to the calculation made using the U-band 
photometry, the dashed curve without it. In all cases 
the photometric redshifts are calculated as the weighted 
average of zPDF as: 

_ J~zP(z)dz 

Zpkot- ^p (z)dz ■ W 

A blank panel in Figure [6] means that the source is out- 
side of the region covered by the CANDELS observations. 

The general result from our analysis is that the photo- 
metric redshifts of the LyC emitter candidates are sys- 
tematically much smaller than the spectroscopic red- 
shifts of their host LBGs. Exceptions include the 
cases of the LyC candidates GDS J033204.91-274451.0, 
CDS J033220.95-275021.8, GDS J033222.95-274727.8, 
and GDS J033223.32-275155.9 (see Figures M and ©, 
where the GOODS photometric redshift is much lower 
than the spectroscopic redshift of the LBG, but the CAN- 
DELS photometric redshift is not. These are cases in 
which the separation between the LBG and the Uem is 
<1" and are not individually detected in the CANDELS 
H-band based catalog (they are marked with red crosses 
in Figure [6]). 

• In the hrst case, GDS J033204.91-274451.0, the LyC 
emitter candidate is resolved as a different source from 
the LBG in the ACS z-band image but not in the H- 
band one. The photometric redshift from the GOODS 
photometry yields a significantly lower redshift than the 
LBG's spectroscopic one, with some dependence of the 
redshift probability function on the exclusion of the U- 
band from the fit, but the resulting photometric redshift 
is in any case much lower than the LBG's. Since the LyC 
candidate and the LBG are not resolved in the H band, 
the CANDELS photometric redshift is essentially that of 
the LBG itself. If the U-band photometry is not included 
in the fit, the spectroscopic and photometric redshift are 
in very good agreement. Including, the U-band data, 
however, yields a photometric redshift somewhat lower 
than the spectroscopic one (z spec = 3.404,z p ^ ot = 3.351), 
which is not surprising, since this is the LBG with the 
lowest redshift in our sample and thus some flux detected 
in the bluest band might easily skew the fit toward lower 
redshift values. 

• The other three sources (GDS J033222.95-274727.8, 
GDS J033220.95-275021.8, GDS J033223.32-275155.9) 
are also cases in which the LyC candidate and the LBG 
are classified as individual sources in the z-band but not 
in the H band. In both cases, the GOODS photometric 
redshift is much lower than the LBG's spectroscopic one. 
In the former, since the LBG redshift is at the high end 
of our sample (z spec = 4.440), adding or not the U-band 



data to the derivation of the photometric redshift makes 
no difference. In the other two, neglecting the U-band 
data results in very good agreement between the spec- 
troscopic and photometric data, while including it lowers 
the photometric redshift solution. 

The LyC emitter candidate with coordinates 
RA=03:32:36.85, DEC=-27:45:57.6, associated with 
galaxy HUDF-J033236.83-274558.0 is discussed in detail 
in the Appendix. It is the fainter among the cases 
discussed in this work with an U mag of 28.6. A new 
observed feature in the J-band (at 10er significance) is 
consistent with relatively strong line emission of a low 
redshift galaxy, lower than its companion LBG. 

A general characteristic of the fits is that including or 
not the U-band photometry in the derivation of the pho- 
tometric redshift of the LyC emitter candidates makes no 
substantial difference in the results and that there is very 
good quantitative agreement between the GOODS and 
CANDELS photometric redshift, pointing to the robust- 
ness of our analysis. 

Finally, a further check has been performed by cal- 
culating the photometric redshifts without the inclusion 
of the U and -B435 bands, that makes the test indepen- 
dent from the single B — V color analysis described in 
Sect. 4.2. Even though the efficiency of the photomet- 
ric redshift prediction is slightly reduced, the resulting 
Zphot best-fit values still favor low redshift solutions for 
all sources, except for GDS J033223. 32-275155.9 that in- 
creases from Zphot = 0.07±0.02 to 3.60±0.2. In this case, 
however, the (B — V) = 0.65 color, the relative escape 
fraction and the A(B — V) measurements are in contrast 
with this high redshift solution. 

4.4. Conclusion From j esc , UV Colors and Photometric 
Redshift Analysis 

In summary, the combination of photometric redshift, 
relative UV colors (that trace the IGM decrement) and 
the relative escape fraction estimated for the LyC can- 
didates all argue against them being ionizing radiation 
escaping form the nearby LBG. Rather, they are most 
likely low-redshift interlopers (as we had originally ar- 
gued in VlOb). The more extreme case among the 
sources here analyzed (HUDF J033236.83-274558.0) is 
discussed in detail in the Appendix, in which we report 
new observational evidence supporting a possible con- 
tamination from a superimposed object. 

Finally, we conclude by noting that deep UV imaging 
with HST/WFC3 will be able to clarify the nature of the 
U-band offset LyC candidates, since the probability that 
LyC radiation at A <C 912 A escapes the host galaxy and 
eludes the IGM is extremely low (for an estimate of this 
probability at A < 700 A see Inoue & Iwata 2008; Meiksin 
2006). 

5. A CRITICAL ANALYSIS OF THE CURRENT 
MEASUREMENTS 

In this section we review all known (to us) galaxies at 
high redshift (z >3) for which emission of LyC ionizing 
radiation has been reported or suspected. 

• Direct detection of LyC ionizing radiation has been 
reported by Shapley et al. (2006, S06) in deep spectra 
of two LBGs at z ~ 3 (U-band dropouts) in the SSA22 
field, dubbed D3 and C49 in their nomenclature, out of a 
total of 14 sources observed with comparable sensitivity. 
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It is worth noting that 109 and Nil found a highly sig- 
nificant null detection for the object SSA22a-D3, which 
is the brighter of the two objects for which S06 claimed 
an ionizing flux detection. 

• Iwata et al. (2009) used deep narrow-band (NB) 
imaging to image the rest-frame SED blueward of the 
912A Lyman limit of members of a well-known overden- 
sity of galaxies at z ~ 3.1, also in the SSA22 field (Steidel 
et al. 1998). They reported 17 detections, 7 from galax- 
ies selected as LBGs and 10 as LAEs. 

• Nestor et al. (2011) used the same technique to ob- 
serve galaxies in the same SSA22 field, but extended the 
sensitivity of previous observations by ss 0.6 mag. They 
reported 34 detections, out of 156 sources, which include 
6 LBGs and 28 LAEs. 

• In the paper VlOb, we used a technique similar, but 
not identical, to the one by 109 and Nil. Instead of using 
a narrow-band, we carried out ultra-deep U-band imag- 
ing to search for candidate LyC emitters among galaxies 
at 3.4 < z < 4.5, namely such that the redshifted 912A 
ionization edge of the targeted sources is to the red and 
completely outside of the filter's bandpass. We found 1 
candidate out of 102 LBGs (described in Sect. 6). 

For most of the NB-selected candidates and LBGs the 
region where the LyC ionization radiation originates is 
found to be spatially offset from that of the non-ionizing 
far-UV light, namely the main body of the galaxy. The 
displacement is generally less than one arcsec, but in 
some cases values as high as 2" have been reported. 
We performed a visual inspection of the images of the 
galaxies by 109 and Nil, and consulted Table 4 and 5 of 
Nil, and found that > 70% of the presently identified 
LyC emitter candidates exhibit a spatial offset between 
the LyC emission and the "non-ionizing" far-UV image 
of the galaxies. Since an offset emission is also consis- 
tent with a lower redshift interloper, a key step is to 
measure the statistical occurrence of these cases (as we 
did in VlOa). In particular with availability of large- 
area HST UV and near-IR surveys with WFC3, such as 
CANDELS, the prospect for substantial progress in this 
regard is good. The other important point is that if the 
offset LyC emission is real, then care must be exercised 
when calculating the fraction of escaping ionizing radia- 
tion / esc in each case, i.e., whether it is from the galaxy 
as a whole or only from a close companion, to avoid bias 
in the derivation of average properties. For example it 
is worth mentioning that the offset possible LyC detec- 
tion in S06, SSA22-C49, has a reported observed ratio 
of (F1500/ F LyC )oBS = 16.4 ±6.1 (Nil), in which the 
flux F1500 used in the calculation is that of the main 
galaxy. Presumably, a calculation restricted to the local 
region where the putative Lyman continuum emission is 
seen would produces a much smaller value, i.e., a much 
higher f esc ,rei- 

Finally, it must be kept in mind that part of the candi- 
dates LyC emitters reported above may have the contri- 
bution to Jose from an AGN component that is difficult 
to identify without multi-wavelength and spectroscopic 
surveys. 

Therefore the amount of stellar LyC emission from 
high-z star-forming galaxies is still uncertain. Very few 
sources have spectroscopic redshift confirmation: (A) 
those (five) reported by Inoue et al. (2011) that need 



extreme stellar populations to justify the flux ratios 
observed, (B) the LBG reported in S06 (SSA22-C49) 
that shows an offset (0.5") U-band emission whose red- 
shift is not conclusively known, and finally (C) the one 
we present in the next section and discovered in the 
GOODS-S field (Ionl). 

6. A LYC EMITTER FOUND IN THE GOODS SAMPLE: 

ION1 

The direct detection of LyC emitters at high redshift 
(as high as possible) , allow us to characterize their prop- 
erties at A > 1215.7A and try to identify similar sources 
and/or look if they are more common during the re- 
ionization epoch. For this reason it is important to iden- 
tify secure LyC emissions. We recall one of the most 
promising stellar ionizers at redshift 3.795 we have iden- 
tified in the GOODS-S field and currently the highest 
redshift known so far (GDS J033216.64-274253.3, Ionl) 
and is useful to revisit it in the context of this work by 
comparing its appearance with the offset Uem sources 
discussed above. 

6.1. Observed Photometric and Spectroscopic Properties 

of Ionl 

This source was initially reported in VlOb, and here 
we briefly summarize and report on new observational 
constraints, and compare it with the offset U em sources 
discussed above. 

SED: The source has not been detected in the F225W, 
F275W and F336W channels of HST/WFG3 observa- 
tions in the GOODS-S (down to 26.3, 26.4, 26.1 at 5a 
for point-like objects, Windhorst et al. (2011)). It lies in 
the B-band dropout selection scheme, with (B-V)=1.70 
and (V-z)=0.40 (V09), and has a blue UV spectral slope, 
having an (i-z) = -0.015, and a p = -2.09 ± 0.16 de- 
rived from Castellano et al. (2011). It is worth noting 
that this source has the bluest (B — V) color among the 
available 17 LBGs with 3.7 < z < 3.9 and compara- 
ble UV slope. This is consistent with the fact that the 
i?435-band is probing the rest-frame interval A < 1020A, 
in which the LyC is contributing to the observed flux 
(a S/N of 10 is measured in this band) and reduces the 
(B — V) color. The probability that zphot > 3.4 is ~ 
100%. It has not been detected in the new 4 Ms Chandra 
X-ray observations, for which we set a la upper limit of 
L x < 3 x 10 42 erg/s, nor in the 24^m SPITZER/MIPS 
observations (Santini et al. 2009). The SED is fully 
comparable to that of a star-forming galaxy, in particu- 
lar from the SED fitting we obtain : SFR ~ 5OM yr _1 , 
stellar mass of 2.3 x 10 9 Af o , E(B - V) < 0.1 (see Fig- 
ure [7]). 

Optical spectrum: Optical spectra from 
Keck/DEIMOS and VLT/VIMOS have been ob- 
tained, although the latter is too shallow and the 
S/N prevent us to add information. From the Keck 
spectrum it is clear that the Lya line is not in emission, 
the continuum-break is evident together with the 
[Si iv] 1393.8-1402.8 and [Civ] 1548. 2-1550.8 absorption 
lines, at redshift 3.795. Also [Cn]1335.1 in absorption 
seems to be present. The other UV absorption lines like 
[Ol] 1302.2, [Si 11] 1260.4 and [Si 11] 1526.7 are not detected 
(see Figure [7]). 

UV and B-band rest-frame morphology: At 
1700A/1900A and B-band rest-frame wavelengths the 
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galaxy shows a resolved and compact shape, with an half- 
light radius of 0.9 kpc (physical size) in both the Zs50 and 
H bands (the latter from CANDELS), respectively (see 
Figure [7|). 

The hypothesis that it is a faint AGN (Lx < 3 x 
10 42 erg/s) is still open. However, the spectral features 
in the optical spectrum and the shape of the SED from 
the U-band to the far infrared (Spitzcr/MIPS 24/xm) are 
fully compatible with a star-forming and relatively low 
mass galaxy. In the following we assume that it is a 
star-forming galaxy. 

The LyC and the nLyC emissions (namely U and 1775- 
bands) are spatially aligned within the errors, A9 < 0.1", 
and the LyC detection at A < 830A implies a rela- 
tively transparent line of sight free from LLSs. From 
Eq. [2] we derive a minimum f esc .rei = 82% assuming 
(L1500/ii a c)int = 7 and a IGM transmission of 0.575 
(maximum value for our U-band filter of the 10000 re- 
alizations at z=3.8). Given its E(B — V) < 0.1 and 
adopting the Calzetti extinction law we have f esc > 56% 
(Calzetti et al. 2000). Assuming a intrinsic value for 
(£1500/ Ll v c) hit = 3 the escape fractions are f eS crei > 
35% and f esc > 24%. 

6.2. How Does Ionl Compare With the Rest of the LBG 
Population ? 

It has been firmly established that as the Lya equiva- 
lent width increases from absorption to emission, the UV 
spectral slope becomes bluer and the strength of the UV 
absorption lines (stellar and interstellar) drastically de- 
creases (see, e.g., Shapley et al. 2003; Kornei et al. 2010; 
Balestra et al. 2010; Vanzella et al. 2009). Moreover, 
on average, high Lya equivalent widths tend to be asso- 
ciated to UV compact morphologies; conversely sources 
with Lya in absorption appears less compact (see Law et 
al. 2007; V09; Pentericci et al. 2010). The Ionl galaxy 
shows : 

1. a blue UV spectral slope (3 = —2.09 ± 0.16, from 
the photometric fitting of Castellano et al. (2011). 

2. relatively weak UV absorption lines, from the 
Keck/DEIMOS spectrum (see Figure [7]). 

3. a compact morphology, 0.9 kpc physical half light 
radius in the 1600A and 3300A rest-frame wave- 
lengths (from HST/ACS and WFC3). 

These are all features that positively correlate with the 
equivalent width of the Lya in emission, nonetheless the 
Lya is absent. In this respect it deviates from the aver- 
age LBG population, particularly if we consider its LyC 
emission. The fact that the Lya line is not in emission 
suggests that it can be absorbed and/or it is intrinsi- 
cally weak because the ionizing radiation is escaping the 
galaxy. The detection of LyC emission would suggest 
low gas attenuation along the line of sight and the blue 
UV slope low dust absorption. In the extreme case in 
which all the ionizing radiation is escaping (/ csc is 100%, 
or is close to), nebular emission lines like Lya , Oxygen 
or the Balmer transitions are no longer pumped, and 
therefore drop to intrinsecally small equivalent widths 
(Schaerer 2003). We can image a geometrical configura- 
tion in which there is a high fraction of escaping ioniz- 



ing radiation and at the same time a high Lya equiva- 
lent width in emission; it can be realized with a sort of 
"unipolar" outflow behind the source that backscatters 
the Lya photons to the observer (as is usually seen), 
while along the line of sight (the front of the galaxy), 
the ISM and CGM are free from dust or gas attenuation 
and the LyC radiation can escapes. However, situations 
like this have never been clearly observed until now and 
the interplay between the escape fraction of ionizing and 
Lya photons is still not known. 

It is worth noting that the absence of Lya emission 
from an object with high / eS c, like Ionl, makes extremely 
difficult the spectroscopic confirmation of similar sources 
at z > 7, i.e. during the reionization epoch. Indeed, 
if place at z = 7 it would be one magnitude fainter, 
Y ~ 25.9, and the continuum break hard to measure 
spectroscopically (Vanzella et al. 2011; Fontana et al. 
2010; Pentericci et al. 2011). 

7. CONCLUSIONS 

The observation of the LyC in distant galaxies is a 
difficult task because several attenuations occur: in the 
galaxy itself and circumgalactic medium (by dust and 
mitral hydrogen gas) and along the intergalactic travel 
(Lyman alpha forest, LLSs and DLAs). A further com- 
plication is related to the possible presence of foreground 
(lower-z) superimposed sources that can mimic the LyC 
emission of the background galaxy. It is therefore crucial 
to establish whether they are genuinely associated with 
the redshift of the main LBG. 

In this work we have discussed three diagnostics that 
can be used to constrain the redshift of these offset and 
faint sources. To this end, the deep and high angular 
resolution multi-wavelength images available from the 
GOODS and CANDELS surveys have been exploited. 
In particular, starting from the sample of U-band emit- 
ters identified in VlOb, we (1) calculated accurate pho- 
tometric redshifts, (2) considered the IGM radial and 
transverse absorption (simulating UV colors and com- 
paring those of LBGs and their U-band emitters) and (3) 
provided new constraints on their nature by calculating 
the f esc quantity under conservative assumption. These 
three analysis suggest that none of offset U-band sources 
are actual sources of LyC emission from the LBGs at 
z > 3.4, but are instead foreground interlopers. This 
strengthens the results of VlOb, i.e., the median f esc 
quantity is very small for L > L* galaxies, or assum- 
ing a bimodal distribution for the f esc , the higher values 
are rare. 

At present, the offset candidates reported in litera- 
ture by S06, 109 and Nil need spectroscopic redshift 
confirmation and are difficult to access if deep multi- 
wavelength high angular resolution imaging is not avail- 
able. 

Moreover, we note that in order to avoid biased mea- 
sures, the calculation of the f esc quantity in the case 
of spatially offset emission must be performed carefully 
considering quantities arising from the same physical re- 
gion. We argue that the measurements reported by S06, 
109 and Nil are affected by this problem and need futhcr 
investigations. 

The current situation is far from clear. The number of 
bona fide LyC detections at z > 3 is very small and not 
statistically significant. 
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We have discovered one good candidate (named Ionl) 
that is currently the highest redshift galaxy known 
with direct LyC detection. Its LyC emission at A < 
830A (probed by the U-band) is aligned with the source 
detected in the UV nLyC, i.e., no offset is observed 
with this resolution. Apart from the three cases re- 
ported by Inoue et al. (2011) with aligned LyC emission 
that need Population-Ill stars to explain their observed 
/ '1500/ / LyC flux ratios (even smaller than one), the Ionl 
emission is easily explainable with standard stellar pop- 
ulation and an f esc > 25%. If the AGN component is 
negligible, as seems to be the case from the X-ray and 
spectroscopic data, it would be the most promising stel- 
lar ionizer that may resemble those responsible for HI 
reionization. 

It is worth noting that a high value of the escape frac- 
tion would correspond to faint nebular emission (e.g., 
Robertson et al. 2010). In this case the Lya line is 
not in emission, despite the fact that the source shows a 
blue UV spectral slope ((3 = —2.09), relatively weak in- 
terstellar absorption lines and compact morphology, all 
characteristics that positively correlate with the strength 
of the Lya emission line (e.g., Shapley et al. 2003; V09; 
Pentericci et al. 2010; Balestra et al. 2010). It is not 
clear if the absence of Lya is connected to the high value 
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of the f esc or due to other effects like the dust attenua- 
tion, that however, would be small, giving the blue UV 
slope. 

Its relatively low stellar mass (2.3 x 10 9 M Q ) may also 
be a signature that feedback processes were more efficient 
to clean the line of sight in this type of galaxy, favoring 
the escape of ionizing radiation. Precisely in this regard 
a measure of the wind velocity is important (Heckman 
et al. 2011; Overzier et al. 2010). We note that Ionl 
is not a LAE, therefore it is important to complement 
the LAE surveys with sources like this. Clearly sources 
like Ionl placed at z>7 would be extremely difficult to 
detect spectroscopically with present facilities. However, 
since at z > 7 the direct measure of LyC is not feasi- 
ble, sources like Ionl are the only viable ways we have 
to investigate the mechanism that allow the ionization 
ratiation to escape, i.e., to address the interplay between 
the f esc quantity and the non-ionizing UV features like 
the strength of the nubular emission lines (e.g., Oxigen, 
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FlG. 1. — Main box: Transmission convolved with the VLT/VIMOS U-band filter (filled circles) as a function of source rcdshift. The 
filled circles and vertical error bars indicate the median value and central 68% range of the transmission for the 10,000 lines of sight 
generated with the new Inoue et al. (2011) simulations. Open squares are the averages calculated over the same lines of sight (shifted by 
dz = 0.05 to the right for clarity). Clearly, the VLT/VIMOS U-band probes progressively shorter wavelengths as rcdshift increases, with 
the effect of lowering the transmission. Inner box: The fraction of recovered sources by the U-band imaging (with S/N > 2) as a function 
of the increasing IGM attenuation (redshift), normalized to the z = 3.4 case (see text for details). 
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Fig. 2. — Examples of 4 LBG (out of 19) with the closest offset U-band detections (and see Figure [8}. If the U— emission arise from a 
region that is at the same redshift as that of the LBG, than it would be genuine LyC emission. The angular separation between the main 
targets (LBG marked in the U-band with 1.2"diameter circles) and the fainter counterparts that generate the U-band emission is smaller 
than 1.0" (see TableQJ. The color ACS/BVI, VIMOS/U and VIMOS/R images are shown for each case, from left to right, respectively, 
and the GOODS ID is reported for the targeted LBG. The box sizes are 6.3" on a side. 
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Fig. 3. — Illustrative example of the measurement of the observed flux ratio between LyC e nLyC for the same system, but in different 
image spatial resolutions. In the top panel the fluxes are measured and compared within the same regions (indicated with green solid 
circles), allowed by the high angular resolution available from HST. In the bottom panels the same is shown for ground-based observations, 
in which the magnitude of the LBG (that includes the fainter offset blob) is wrongly adopted as /1500, as it has been done in recent works. 
The two ratios differ by more than a factor of 10 (see text and Table [l] for details) like their / csc rc i. The solid green circles of the top 
panels are reported also in the bottom-left (dashed circle) and right panels. The black solid circles indicate the position of the LBG. In 
this example we consider negligible the difference between the nLyC fluxes probed by the 1775 and R bands. 
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Fig. 4. — Differences between the UV colors of LBGs and the U-band emitters (A(_B — V) = (B — V)lbg — (B — V)jjem) as a function 
of redshift (bottom) and physical transverse separation of the pairs (LBG-Uemi at the mean redshift of the Lyman-alpha forest (top). 
Open circles are the observed differences, whose typical error is < 0.2 (see Tab. [l] and their S/N in the £?435-band). Solid circles are the 
same quantity after correcting the colors by running MC simulations: a Lyman continuum emission with f esc = 1 has been inserted in the 
£?435-band of the Uem. Filled triangles are the colors derived from the same MC simulations by assuming a conservative IGM transmission 
of 100% (see text for details). Green symbols represents the A(V — i) for the pairs with LBG at 4.0 < z < 4.5. In these cases the V606 _ t )an d 
is more suitable than the i?435-band and includes only the IGM at A > 912A (no LyC is included in the V606 filter (see text)). 
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Fig. 5. — Photometric redshift probability function of the LyC emitter candidates derived from the GOODS photometry in the 
VLT/VIMOS U, HST/ACS BViz, VLT/ISAAC JHK and Spitzer/IRAC h, h, h, U bands in matched apertures measured with the 
TFIT photometric package (see text). The dashed curve shows the case where the U band photometry is not used in the calculation. The 
galaxy with the blank panel (and Z p (j Ot =-1.0) is the case described in Appendix. The two Zp^ot values (top and bottom) in each panel 
indicate the redshift calculated with and without the U-band, respectively. The Z^bq reports the spectroscopic redshift (top) and the 
photometric redshift (bottom). In the bottom right panel the Ionl source with LyC detection is shown. 
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FlG. 6. — Photometric redshift probability function of the LyC emitter candidates derived from the CANDELS photometry in the 
VLT/VIMOS U, HST/ACS BViz, HST/WFC3 YJH and Spitzer/IRAC 7i, 7 2 , I3, ^4 bands in matched apertures measured with the TFIT 
photometric package (see text). The dashed curve shows the case where the U band photometry is not used in the calculation. A galaxy 
with a blank panel (and Z p h ot =-l.Q) means that no WFC3 data are available for that source. The two Z^^t values (top and bottom) 
ineach panel indicate theredshift calculated with and without the U-band, respectively. The Z^bg reports the spectroscopic redshift (top) 
and the photometric redshift (bottom). In the bottom right panel the lonl source with LyC detection is shown. Crosses are the Uem 
sources for which the photometric redshift is not useful (see text for further details). 
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lonl : J033216.64-274253.3, redshift 3.795 




wavelength (A) X-ray 4Ms VIMOS/U-band HSTMCS-BVi HST/WFC3 

Chandra >.<830A GOODS F160W 

CAN DELS 



Fig. 7. — SED fitting (left, panel A), Keck/DEIMOS spectrum (top right, panel B) and image cutouts of the lonl galaxy (bottom right, 
panel C) are shown. Panel A: The SED fitting is shown and includes from left to right the VLT/VIMOS [/-band, HST/ ACS B435, V 6 06, 
*775, 2850 bands and the HST/WFC3 Y, J and H bands, the VLT/ISAAC A"g-band (the most deviant point), and the four Spitzer /IRAC 
channels 3.6fim, 4.5/wra 5.8^tm and 8.0/wra (only two points come from ground-based observations). Panel B: Keck/DEIMOS spectrum 
with UV absorption features indicated that yield a redshift z = 3.795 (solid vertical lines and labels). Dotted vertical lines mark the 
expected positions of [Ol]1302.2, [Sin]1260.4 and [Si II] 1526.7, from left to right, respectively. Panel C: Image cutouts from the X-ray to 
the U— band rest-frame are shown (the VIMOS, ACS and WFC3 box sizes and the square inset in the X-ray image are 6.3" on a side). 
Note that the VLT/VIMOS [/-band is sampling the wavelengths smaller than 830A rest-frame. 
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TABLE 1 

Sample of offset U-band emission from LBG in GOODS-S. 
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/ /1500 \ 
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f a 

•> esc, ret 

(%) 
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1.20 
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1.98(0.23),-1.00 


19,16 


1.71 
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0.66 


792 
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1.50 


4.292 


1.44(0.09), 1.48(0. 05) 


30, <1 
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1.15 


1.01 


724 


J033217. 47-274141. 6 


1.20 


3.780 


0.52(0.04),0.51(0.03) 


20,4 


1.00 


— 


6.85 


76 


J033217. 41-275200. 3 


1.40 


3.494 


0.93(0. 18), 0.67(0. 24) 


10,15 


1.07 


- 


2.29 


197 


J033217. 79-275050. 4 


1.60 


3.477 


1.55(0. 03), 1.52(0. 08) 


8,13 
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TA0001 o or o^r 1 o a o 

J033218. 85-275134. 3 
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3.659 
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2.17 
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o. 1 04 


i on/n 1 a \ i 7K(r\ m\ 
l.yu^U. 14 J, 1. 1 o^U.Uo j 


on ^ 1 


A.oA 




1.04 


ooy 


J033220.95-275021.8* 


0.60 
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2.66(0.07),-1.00 
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3.08 


116 


J033222. 79-274439.0 


1.30 


3.582 


0.58(0.75),0. 39(0.02) 


28,11 


0.85 




2.27 


210 


J033222. 95-274727.8* 


0.80 


4.440 


0.46(0.46),-1.00 


4,<1 




1.05 


2.31 


325 


J033223. 32-275155.9 


1.10 


3.470 


0.07(0.07),-1.00 
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1.36 




4.45 
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J033226. 43-274124.5* 


0.70 


4.384 


2.31(0.16), 2.32(0.20) 


19,4 
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J033228. 95-274232. 9 


1.40 


3.585 


2.72(0.18), 2. 77(0. 17) 


12,11 


1.05 




2.42 


197 


J033234.33-274123.0 


1.60 


3.418 


2.30(0.25), 2. 33(0. 31) 
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J033236. 85-274557.6' 
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1.36 


387 


J033240.48-274431.7 


1.40 


4.120 


2.39(0.07), 1.87(0.04) 
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0.54 
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1.10 


3.472 


2.17(0.25),-1.00 
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1.06 




1.66 


272 


(Ionl) 


<0.1" 


3.795 


3. 64(0. 04), 3. 68(0. 09) 


10.3 


15.00 


35 



Note. — (a) These are the minimum f eac re j derived from Eq. [2] adopting conservatively the maximum IGM transmission at the given 
redshift. 10000 IGM transmissions have been calculated from simulations specifically computed at the spectroscopic redshift of the LBG. 
The intrinsic ratio has been assumed to be L1500/ Ll v c = 3 (see text for details), (b) galaxy observed in the HUDF and shown in Figure|8] 
The angular separation between the compact blue emitter and the center of light of the LBG is 0.4". (c) The photometric redshifts with 
their r.m.s. (GOODS and CANDELS, on left and right, respectively) without the inclusion of the U-band are reported. (*) sources shown 
in Figure [2] 



APPENDIX 

THE CASE OF THE LBG HUDF .1033236.83-274558.0 

Figure |S] shows the HST/ ACS images of the offset U-band LyC emitter candidate that we have identified at 0.4" 
from LBG HUDF J033236.83-274558.0 at z = 3.797 (its celestial coordinates are RA=03:32:36.85, DEC=-27:45:57.6). 
This pair has already been discussed by VlOa, who argued that the LyC candidate is most likely a low redshift 
interloper. Here we add another piece of evidence in support of this interpretation. We have complemented the ACS 
optical photometry, which takes advantage of the HUDF images (Bcckwith et al. 2006), with the new WFC3/IR 
data obtained as part of the CANDELS project (Koekemoer et al. 2011; Grogin et al. 2011) in the Y (F105W), 
J (F125W) and H (F160W) bands. The SED of the Uem from 0.5 to 1.6 \im is blue, essentially consistent with a 
flat spectrum, i.e., /„ oc v , except in the J band where the source is ~ 0.43 mags brighter (at 10-sigma), namely 
C/=28.63±0.2, J B 435 = 28.59±0.05, U 606 =28.31±0.03, i 775 =28.31±0.03 , 2 850 = 28.23 ±0.05, F=28.31±0.07, J=27.88±0.04 
and iJ=28.32±0.06. Using the HST photometry alone, with or without the U band, the photometric redshift remains 
unchanged, i.e., z p hot — 3.8. 

At the redshift of the LBG, z — 3.797, the J band probes the rest-frame spectral range AA - 2295, 2918 A where 
there are no strong emission lines, except for the possible [Mgii] 2797,2803 P-Cygni feature, which, however, has not 
been observed to yield any net emission in star-forming galaxies at redshift up to z ~ 2 (e.g., Weiner et al. 2009; 
Rubin et al. 2010; Giavalisco et al. 2011). The simplest and, in our view, most convincing interpretation is that the 
U-band source is at lower redshift and the J band flux is boosted by an emission line (or more than one), most likely 
[On] or H/3 and [Oiii]. 

All photometric redshift solutions lie at z < 3. The required rest-frame equivalent width of a single line to boost 
the J-band flux of 0.4 magnitude is 670A if Ha at z ~ 0.9 (L Ha = 1.5 x 10 40 er 5 /s), or 383A if [On]3727 at z ~ 2.5 
(L [OII] = 1.9 x 10 41 er 5 / S ). 

Also the complex of the three lines H/3, [O III] 4959 and [O in] 5007 can fall in the suitable wavelength interval. In 
the case of Ha we would detect [On]3727, H/3, and [O in] 4959 - 5007 in the spectrum of the LBG, that is not the 
case. Similarly, in the case of H/3 and [O in] complex, the [O n]3727 line would be present in the spectrum of the LBG 
at ~ 9700A. However in this latter case it may be too faint to be detected. Therefore we exclude the Ha possibility 
only, and keep the other two configurations. 

Given the response curve of the WFC3/Y/J/H filters, in order to have the above complex of three lines or the single 
[O n]3727 line in the J filter only (green shaded region of Figure [8|), the redshift of the source would be 1.45 < z < 1.8 
or 2.25 < z < 2.7, respectively. 

As mentioned above, the required rest-frame equivalent widths are high (larger than 300A), however we note that 
ultra-strong emission line (compact) galaxies with these features have been recently detected (Kakazu et al. 2007; Hu 
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et al. 2009; Izotov et al. 2011 and more recently van der Wei et al. 2011). Our source would extend their findings to a 
more extreme limit of luminosity, i.e., a rest-frame absolute i?435-band magnitude of -16.86. Going a bit more into the 
details, we note that the SED shows a significant decrease of flux in the U and -B435 bands (see Figure [8]), that would 
not be consistent with the redshift z ~ 1.6, which corresponds to A > 1200A. Conversely, in the case of [Oll]3727 at 
z ~ 2.5, these bands are affected by the IGM and the Lyman limit attenuation. Indeed the best photometric redshift 
solution is zphot ~ 2.5. Therefore we argue that the compact source describe here is at z < 3 and a probable strong 
Oxygen emitter (similar to those recently identify in the CANDELS survey by van der Wei et al. 2011), that mimic a 
LyC emission of the background LBG. 
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Fig. 8. — Images and spectral energy distribution for J033236. 83-274558.0 in the HUDF. In the top part, from left to right: the deep 
color composite obtained with H5T/WFC3 F105 {Y), F125 (J), F160 (H) bands, the BVI color HST/ACS composite, the ultra-deep 
VLT/VIMOS U-band image (that is probing the LyC for the LBG) and the deep VLT/VIMOS R band image. The box sizes are 4.5" on a 
side. Note that such a system observed from the ground is poorly resolved (ij-band on the right). The dotted arrows indicate the position 
of the LBG at z = 3.797 and the solid arrows mark the position of the blue emitter. The angular separation between the blue emitter and 
the peak of the LBG emission is ~ 0.4". However the LBG extends beyond this separation and the blue spot is fully superimposed. This 
compact blue source shows also a boosted emission in the J band, ~ 0.4 magnitudes brighter (10<r) than the other two Y and H (bottom 
left), it is the reason why it appears as a green spot in the WFC3 composite. This J-band emission is compatible with a low redshift 
interpretation (zphot < 3, see text for details). In the bottom right, the ffi>T/WFC3 relative transmissions of the Y, J and H niters are 
shown. The shaded green region mark the wavelength range where the presence of the emission line(s) are boosting the J-band flux. 



